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Abstract 

The longstanding open problem of approximating all singular vertex couplings in a quantum 
graph is solved. We present a construction in which the edges are decoupled; an each pair of 
their endpoints is joined by an edge carrying a 6 potential and a vector potential coupled to the 
"loose" edges by a (5 coupling. It is shown that if the lengths of the connecting edges shrink 
to zero and the potentials are properly scaled, the limit can yield any prescribed singular vertex 
coupling, and moreover, that such an approximation converges in the norm-resolvent sense. 
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1. Introduction 

While the origin of the idea to investigate quantum mechanics of particles confined to a graph 
was conceived originally to address to a particular physical problem, namely the spectra of aro- 
matic hydrocarbons 1 1 1, the motivation was quickly lost and for a long time the problem remained 
rather an obscure textbook example. This changed in the last two decades when the progress of 
microfabrication techniques made graph-shaped structures of submicron sizes technologically 
important. This generated an intense interest to investigation of quantum graph models which 
went beyond the needs of practical applications, since these models proved to be an excellent 
laboratory to study various properties of quantum systems. The literature on quantum graphs is 
nowadays huge; we limit ourselves to mentioning the recent volume |2| where many concepts 
are discussed and a rich bibliography can be found. 

The essential component of quantum graph models is the wavefunction coupling in the ver- 
tices. While often the most simple matching conditions (dubbed free, Kirchhoff, or Neumann) 
or the slightly more general S coupling in which the functions are continuous in the vertex are 
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used, these cases represent just a tiny subset of all admissible couplings. The family of the latter 
is determined by the requirement that the corresponding Hamiltonian is a self-adjoint operator, 
or in physical language, that the probability current is conserved at the vertices. It is not difficult 
to find all the admissible conditions mathematically; if the vertex joins n edges they contain 
free parameters, and with exception of the one-parameter subfamily mentioned above they are 
all singular in the sense that the wavefunctions are discontinuous at the vertex. 

What is much less clear is the physical meaning of such conditions. It is longstanding open 
problem whether and in what sense one can approximate all the singular couplings by regular 
ones depending on suitable parameters, and the aim of the present paper is to answer this question 
by presenting such a construction, minimal in a natural sense using real parameters, and to 
show that the closeness is achieved in the norm-resolvent sense, so the convergence of all types 
of the spectra and the corresponding eigenprojections is guaranteed. 

The key idea comes from a paper of one of us with Shigehara [31 which showed that a com- 
bination of regular point interactions on a line approaching each other with the coupling scaled 
in a particular way w.rt. the interaction distance can produce a singular point interaction. Later 
it was demonstrated [4] that the convergence in this model is norm-resolvent and the scaling 
choice is highly non-generic. The idea was applied by two of us to the simplest singular cou- 
pling, the so-called 6'^, in |5| and was demonstrated to work; the question was how much it can 
be extended. Two other of us examined it |6| and found that with a larger number of regular 
interactions one can deal with families described by 2n parameters, and changing locally the 
approximating graph topology one can deal with all the couplings invariant with respect to the 
time reversal which form an ("2')-parameter subset. 

It was clear that to proceed beyond the time-reversal symmetry one has to involve vector 
potentials similarly as it is was done in the simplest situation in Q. In this paper we present 
such a construction which contains parameters breaking the symmetry and which at the same 
time is more elegant than that of [6J in the sense that the needed "ornamentation" of the graph 
is minimal: we disconnect the n edges at the vertex and join each pair of the so obtained free 
ends by an additional edge which shrinks to a point in the limit. The number of parameters leans 
on the decomposition rp- - n + 2(2), where the first summand, «, corresponds to 6 couplings of 
the "outer" edge endpoints with those of the added shrinking ones. The second summand can 
be considered as times two parameters: one is a 5 potential placed at the edge, the other is a 
vector potential supported by it. 

Our result shows that any singular vertex coupling can be approximated by a graph in which 
the vertex is replaced by a local graph structure in combination with local regular interactions and 
local magnetic fields. This opens way to constructing "structured" vertices tailored to the desired 
conductivity properties, even tunable ones, if the interactions are controlled by gate electrodes, 
however, we are not going to elaborate such approximations further in this paper. 

We have to note for completeness that the problem of understanding vertex couplings has also 
other aspects. The approximating object needs not to be a graph but can be another geometrical 
structure. A lot of attention was paid to the situation of "fat graphs", or networks of this tubes 
built around the graph skeleton. The two approaches can be combined, for instance, by "lifting" 
the graph results to fat graphs. In this way approximations to 5 and 6'^ couplings by suitable 
families of Schrodinger operators on such manifolds with Neumann boundaries were recently 
demonstrated in |8|. The results of this paper can be similarly "lifted" to manifolds; that will be 
the subject of a subsequent work. 

Let us review briefly the contents of the paper. In the next section we gather the needed 
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preliminary information. We review the information about vertex couplings and derive a new 
parametrization of a general coupling suitable for our purposes. In Section 3 we describe in 
detail the approximation sketched briefly above and show that on a heuristic level it converges 
to a chosen vertex coupling. Finally, in the last section we present and prove our main result 
showing that the said convergence is not only formal but it is valid also in the norm-resolvent 
sense. 



2. Vertex coupling in quantum graphs 

Let us first recall briefly a few basic notions; for a more detailed discussion we refer to the 
literature given in the introduction. The object of our interest are Schrodinger operators on metric 
graphs. A graph is conventionally identified with a family of vertices and edges; it is metric if 
each edge can be equipped with a distance, i.e. to be identified with a finite or semi-infinite 
interval. 

We regard such a graph F with edges £1 ,...,£■„ as a configuration space of a quantum me- 
chanical system, i.e. we identify the orthogonal sum = 0"^j L\Ej) with the state Hilbert 
space and the wave function of a spinless particle "living" on F can be written as the column 
^ - (i/'i , \l/2, • ■ ■ , ijJnY with ijjj e L?-{E j). In the simplest case when no external fields are present 
the system Hamiltonian acts as {HY^)j - -^'p with the domain consisting of functions from 
W^''^{T) := 0'L, W^'^iEj). Not afl such functions are admissible, though, in order to make the 
operator self-adjoint we have to require that appropriate boundary conditions are satisfied at the 
vertices of the graph. 

We restrict our attention to the physically most interesting case when the boundary conditions 
are local, coupling values of the functions and derivatives is each vertex separately. Our aim 
is explain the meaning of a general vertex coupling using suitable approximations; the local 
character means that we can investigate how such a system behaves in the vicinity of a single 
vertex. A prototypical example of this situation is a star graph with one vertex in which a finite 
number of semi-infinite edges meet; this is the case we will mostly have in mind in the following. 

Let us thus consider a graph vertex V of degree n, i.e. with n edges connected at V. We 
denote these edges by Ei, . . . ,E„ and the components of the wave function values at them by 
^i{xi), . . . , ij/„{x„). We choose the coordinates at the edges in such a way that xj > for all 
j - and the value xj = corresponds to the vertex V. For notational simplicity we put 

= (lAi (0), ■ ■ ■ , 4'n{'^)Y ™d '^'y - (i/r'j (0), . . . , 4>',SS^)Y ■ Since our Hamiltonian is a second-order 
differential operator, the sought boundary conditions will couple the above boundary values, their 
most general form being 

A^y + BYy^O, (1) 

where A and B are complex n x n matrices. 

To ensure self-adjointness of the Hamiltonian, which is in physical terms equivalent to con- 
servation of the probability current at the vertex V, the matrices A and B cannot be arbitrary but 
have to satisfy the following two conditions, 

• rank(A|B) = n, 

, (2) 

• the matrix AB is self-adjoint, 

where {A\B) denotes the n x 2n matrix with A,B forming the first and the second n columns, 
respectively, as stated for the first time by Kostrykin and Schrader [9J. The relation ([1]) together 
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with conditions (|2]l (for brevity, we will write ([T]i&(|2]l) describe all possible vertex boundary 
conditions giving rise to a self-adjoint Hamiltonian; we will speak about admissible boundary 
conditions. 

On the other hand, it is obvious that the formulation ([T]i&(|2]) is non-unique in the sense that 
different pairs (Ai, Bi), (A2, B2) may define the same vertex coupling, as A, B can be equivalently 
replaced by CA,CB for any regular matrix C e C"". To overcome this ambiguity, Harmer 
HP I, and independently Kostrykin and Schrader IITTI proposed a unique form of the boundary 
conditions ([T]), namely 

(U-I)^'v + S,(U + IWv^O, (3) 

where f/ is a unitary n x n matrix. Note that in a more general context such conditions were 
known before (TT\, see also (13]. 

The natural parametrization ([3]) of the family of vertex couplings has several advantages in 
comparison to ([TJ&Q, besides its uniqueness it also makes obvious how "large" the family is: 
since the unitary group U{n) has n^ real parameters, the same is true for vertex couplings in 
a quantum graph vertex of the degree n. Of course, this fact is also clear if one interprets the 
couplings from the viewpoint of self-adjoint extensions lfT4l . 

On the other hand, among the disadvantages of the formulation ([3]) one can mention its com- 
plexity: vertex couplings that are simple from the physical point of view may have a complicated 
description when expressed in terms of the condition (|3]l. As an example, let us mention in the 
first place the 5-coupling with a parameter a e M, characterized by relations 



(A/O) = iA*(0) «A(0) , =!...,«, ^ .A>(0) = t^^^(O) , (4) 
for which the matrix U used in (|3]l has entries given by 



Ujk = ^ - Sjk , (5) 

5jk being the Kronecker delta. When we substitute (|5]) into (|3]l and compare with Q rewritten 
into a matrix form ([TJ, we observe that the first formulation is not only more complicated with 
respect to the latter, but also contains complex values whereas the latter does not. This is a reason 
why it is often better to work with simpler expressions of the type ([T])&(|2|. Another aspect of this 
parametrization difference concerns the meaning of the parameters. Since the n^ ones mentioned 
earlier are "encapsulated" in a unitary matrix, it is difficult to understand which role each of them 
plays. 

On the other hand, both formulations ([T]i&(|2]l and Q have a common feature, namely that 
they have a form insensitive to a particular edge numbering. If the edges are permuted one has 
just to replace the matrices A, B and U hy A,B and U, respectively, obtained by the appropriate 
rearrangement of rows and columns. This may hide different ways in which the edges are cou- 
pled; it is easy to see that a particular attention should be paid to "singular" situations when the 
matrix U has eigenvalue(s) equal to +1. 

Since the type of the coupling will be important for the approximation we are going to con- 
struct, we will rewrite the vertex coupling conditions in another form which is again simple and 



unique but requires an appropriate edge numbering. This will be done in Theorem 2.1 before 
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stating it we introduce several symbols that will be employed in the further text, namely 



C*'' - the set of complex matrices with k rows and / columns, 

h - the set {1,2, . . . ,n), 
/(«) _ the identity matrix n X n. 

To be precise, let us remark that the term "numbering" with respect to the edges connected in the 
graph vertex of the degree n means strictly numbering by the elements of the set «. 

Theorem 2.1. Let us consider a quantum graph vertex V of the degree n. 

(i) Ifm < n, S € C'"'"' is a self-adjoint matrix and T e C'"'""™, then the equation 



j(m) 





S 

j(n-m) 



(6) 



expresses admissible boundary conditions. This statement holds true for any numbering 
of the edges. 

( ii) For any vertex coupling there exist a number m < n and a numbering of edges such that 
the coupling is described by the boundary conditions (|6]l with the uniquely given matrices 
T e C"''"-'" and self-adjoint S e 

( Hi) Consider a quantum graph vertex of the degree n with the numbering of the edges explicitly 
given; then there is a permutation 11 G 5„ such that the boundary conditions may be written 
in the modified form 



j{m) 





s 





j(n-m) 



(7) 



for 





' <^n(i)(0) ^ 




' 'An(,)(0) ^ 












. iAn(«)(0) , 




. '/'n(„)(0) , 



where the self-adjoint matrix S e C"'™ and the matrix T e C"*'""™ depend unambigu- 
ously on n. This formulation of boundary conditions is in general not unique, since there 
may be different admissible permutations 11, but one can make it unique by choosing the 
lexicographically smallest permutation 11. 

Proof. The claim (iii) is an immediate consequence of (ii) using a simultaneous permutation of 
elements in the vectors and *l"y, so we have to prove the first two. As for (i), we have to show 
that the vertex coupling ([T]) with matrices 



S 

T* 





_j-{n-m) 



and B 



conform with (|2|. We have 
rank 



-S 



/C™) T 

j(n-m) Q Q 



rank 
5 



/(«') 




/(«) Y 






_j{n-m) 



-S T 

T* 



and 

-S W /<«' T \* I -S \ 
T* _/("-'«) )'[ j ~\ oj' 

the latter matrix is self-adjoint since S - S*, thus (|2]i is satisfied. 

Now we proceed to (ii). Consider a quantum graph vertex of the degree n with an arbitrary 
fixed vertex coupling. Let and ^''y denote the vectors of values and derivatives of the wave 
function components at the edge ends; the order of the components is arbitrary but fixed and the 
same for both vectors. We know that the coupling can be described by boundary conditions ([T]l 
with some A,B e C"'" satisfying (|2]). Our aim is to find a number m < n, a certain numbering 
of the edges and matrices S and T such that the boundary conditions ([T]i are equivalent to (|6|. 
Moreover, we have to show that such a number m is the only possible and that S,T depend 
uniquely on the edge numbering. 

When proceeding from ([T]) to (|6|, we may use exclusively manipulations that do not afifect 
the meaning of the coupling, namely 

• simultaneous permutation of columns of the matrices A, B combined with corresponding 
simultaneous permutation of components in ^Py and 

• multiplying the system from left by a regular matrix. 

We see from (|6]) that m is equal to the rank of the matrix applied at We observe that the 
rank of this matrix, as well as of that applied at *Fv, is not influences by any of the manipulations 
mentioned above, hence it is obvious that m - rank(B) and that such a choice is the only possible, 
i.e. m is unique. 

Since rank(B) = m with me {0, . . . , n), there is an m-tuple of linearly independent columns of 
the matrix B; suppose that their indices are . . , jm- We permute simultaneously the columns 
of B and A so that those with indices ji, . . . , jm are now at the positions 1, . . . , m, and the same 
we do with the components of the vectors ^'y. Labelling the permuted matrices A, B and 
vectors ^'y with tildes, we get 

A^'v + B% = . (8) 

Since rank(B) = rank(B) - m, there are m rows of B that are linearly independent, let their 
indices be /i, . . . , and n - m rows that are linear combinations of the preceding ones. First we 
permute the rows in (|8]l so that those with indices ii, . . . , i,„ are put to the positions 1, . . . , m; note 
that it corresponds to a matrix multiplication of the whole system ([8]l by a permutation matrix 
(which is regular) from the left, i.e. an authorized manipulation. In this way we pass from A 
and B to matrices which we denote as A and B; it is obvious that this operation keeps the first m 
columns of the matrix B linearly independent. 

In the next step we add to each of the last n - m rows of AV(Q) + B^'{Q) = such a linear 
combination of the first m rows that all the last n - m rows of B vanish. This is possible, because 
the last n - m lines of B are linearly dependent on the first m lines. It is easy to see that it 
is an authorized operation, not changing the meaning of the boundary conditions; the resulting 
matrices at the LHS will be denoted as B and A, i.e. 

A^'v + B% = . (9) 

From the construction described above we know that the matrix B has a block form, 

j' 
6 



where Sn e C™'™ and S12 e ([yhn-m.^ square matrix Sn e C'"'"' is regular, because its 
columns are linearly independent. We proceed by multiplying the system (|9| from the left by the 
matrix 



-I 
II 

Q j(n-m) 



'u 



arriving at boundary conditions 

^1 

1^21 ■^ii I \ 
where S12 = S7/Si2. 

Boundary conditions ( [TO] i are equivalent to ([T]l, therefore they have to be admissible. In other 

I and „ have to satisfy both the conditions (|2|l, 

-^21 -^22 / \ ^ ^ / 
which we are now going to verify. Let us begin with the second one. We have 

/ ^11 \ ( \ / ^„ +Jli2S\2 \ 

\ ^21 ^22 j \ ^21 +^22^*2 j 

and this matrix is self-adjoint if and only if + ^12^*12 ^^^^ adjoint and ^2\ + ^22^X2 - 0- 
We infer that ^(21 = -^22^12^ hence condition ([TO]l acquires the form 



^11 ^12 , / 



-^^22^;, y.22 j^^'M'o 0" = ^11) 



The first one of the conditions (|2]i says that 



yill yii2 



hence xwk[-m22'B\^\^22) = n-m. Since (-yi22S';2lJ?l22) = -^22 ■ we obtain the 

condition rank(J?l22) - n — m, i.e. ^22 must be a regular matrix. It allows us to multiply the 
equation ( [TT| from the left by the matrix 

/C") -^12^22 

~^22 

which is obviously well-defined and regular; this operation leads to the condition 

-/("-"') j^^' + l^ j^^ 

If follows from our previous considerations that the square matrix + J^nSlj is self-adjoint. 
If we denote it as -S , rename the block S12 as T and transfer the term containing "^'y to the right 
hand side, we arrive at boundary conditions 

T \,.., I S 



^ \ -T* /(" 



"Py - T-t jln-m) ■ (12) 
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The order of components in '¥v and "^'y determines just the appropriate numbering, in other 
words, the vectors ^'y and 'i'y represent exactly what we understood by ^'y and ^''y in the for- 
mulation of the theorem. 

Finally, the uniqueness of the matrices S and T with respect to the choice of the permutation 
n is a consequence of the presence of the blocks /^'"^ and First of all, the block '"^ 

implies that there is only one possible T, otherwise the conditions for i/','„_|_i , ■ ■ ■ ,4''n would change, 
and next, the block 7*™' together with the uniqueness of T implies that there is only one possible 
S , otherwise the conditions for ipi, . . . , {}/,„ would change. □ 

Remark 2.2. The expression (|7]i implies, in particular, that if B has not full rank, the number of 
real numbers parametrizing the vertex coupling ([T]) is reduced from rp- to at most m{2n - m) - 
rp- - in - nif-, where m - rank(B). Another reduction can come from a lower rank of the matrix 
A. 

Remark 2.3. The procedure of permuting columns and applying linear transformations to the 
rows of the system ([!]) has been done with respect to the matrix B, but one can start by same 
right from the matrix A as well. In this way we would obtain similar boundary conditions as (|6|, 



only the vectors ^'y and T'^ would be interchanged. Theorem 2.1 can be thus formulated with 
Equation (|6| replaced by 

//W T\ I S \ 
\ / \ -T* /(«-'") I V ■ 

For completeness' sake we add that another possible forms of Equation (|6]l in Theorem 2. 1 

{ S \ / /('") T\ 
I .T'* /("-'«) 1^1 I ^ ~ 

and 

r \,„ Is 



-T* /'«-'") j^'v"*^' 

having the standardized form A^y + B'V'y - 0, last two formulations may be sometimes more 
convenient than (|6|. 

Obviously, an analogous remark applies to Equation (|7]i. 

Remark 2.4. A formulation of boundary conditions with a matrix structure singling out the 
regular part as in (|7]i has been derived in a different way by P. Kuchment |15|. Recall that in 
the setting analogous to ours he stated existence of an orthogonal projector P in C" with the 
complementary projector Q - Id-P and a self-adjoint operator L in Q€." such that the boundary 
conditions may be written in the form 

"^^^ ^ ' (13) 

Let us briefly explain how P. Kuchment's form differs from (|7]). When transformed into a matrix 
form, ([OJ consists of two groups of n linearly dependent equations. If we then naturally extract 
a single group of n linearly indepent ones, we arrive at a condition with a structure similar to 
( [TT| , i. e. the upper right submatrix standing at ^''y is generally a nonzero matrix mx (n - m). 
In other words, whilst P. Kuchment aimed to decompose the boundary conditions with respect to 



two complementary orthogonal projectors, our aim was to obtain a unique matrix form with as 
many vanishing terms as possible; the form (|6| turned out to have a highly suitable structure for 
solving the problem of approximations that we are going to analyze in the rest of the paper 

To conclude this introductory section, let us summarize main advantages and disadvantages 
of the conditions (|6]) and Q. They are unique and exhibit a simple and clear correspondence 
between the parameters of the coupling and the entries of matrices in (|6]l, furthermore, the matri- 
ces in (|6]l are relatively sparse. On the negative side, the structure of matrices in (|6]l depends on 
rank(B) and the vertex numbering is not fully permutable. 

3. The approximation arrangement 

We have argued above that due to a local character one can consider a single-vertex situation, 
i.e. star graph, when asking about the meaning of the vertex coupling. In this section we consider 
such a quantum graph with general matching conditions and show that the singular coupling may 
be understood as a limit case of certain family of graphs constructed only from edges connected 
by ^-couplings, (5-interactions, and supporting constant vector potentials. 

Following the above discussion, one may consider the boundary conditions of the form (|6|, 
renaming the edges if necessary. It turns out that for notational purposes it is advantageous to 
adopt the following convention on a shift of the column indices of T: 

Convention 3.1. The lines of the matrix T are indexed from 1 to m, the columns are indexed 
from m + \ to n. 

Now we can proceed to the description of our approximating model. Consider a star graph 
with n outgoing edges coupled in a general way given by the condition (|7]). The approximation 
in question looks as follows (cf. Fig.l): 

• We take n halflines, each parametrized by jc e [0, H-oo), with the endpoints denoted as Vj, 
and put a ^-coupling (to the edges specified below) with the parameter Vj{d) at the point 
Vj for all yen. 

• Certain pairs Vj, Vk of halfline endpoints will be joined by edges of the length 2d, and 
the center of each such joining segment will be denoted as W\j,k]- For each pair {j,k], the 
points Vj and Vk, j + k, are joined if one of the following three conditions is satisfied (keep 
in mind Convention |3.1[ ): 

(1) j e m, k > m + I, and Tj^ (or j > m + I, k e m, and T^j + 0), 

(2) j,k€m and {31 > m + \){Tji i^OATui^ 0), 

(3) j, k € m, S jii + 0, and the previous condition is not satisfied. 

• At each point we place a d interaction with a parameter 'W\j}i\{d). From now on we 
use the following convention: the connecting edges of the length 2d are considered as 
composed of two line segments of the length d, on each of them the variable runs from 
(corresponding to the point ^^[jji]) to d (corresponding to the point V j or V^). 

• On each connecting segment described above we put a vector potential which is constant 
on the whole line between the points V, and V^- We denote the potential strength between 
the points and V j as A(jk)(d), and between the points VK|j,jt) and V^ as A(itj)(d). It 
follows from the continuity that A^^jyid) = -A^j^^iid) for any pair {j, k). 
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Figure 1: The scheme of the approximation. All inner links are of length 2d. Some connection links may be missing if 
the conditions given in the text are not satisfied. The quantities corresponding to the index pair (j, k} are marked, and the 
grey line symbolizes the vector potential A^j^|^^{d). 



The choice of the dependence of Vj{d), W{j^k){d) and A(j;jt)(<i) on the parameter d is crucial for the 
approximation and will be specified later. 

It is useful to introduce the set Nj c h containing indices of all the edges that are joined to 
the 7-th one by a connecting segment, i.e. 



For the wave function components on the edges we use the following symbols: 

• the wave function on the 7-th half Une is denoted by if/j, 

• the wave function on the line connecting points Vj and Vk has two components: the one on 
the line between W^j^^ and Vj is denoted by <P{ik)^ the one on the half between the middle 
and the endpoint of the A:-th half Une is denoted by ip{k,j)- We remind once more the way in 
which the variable x of ip(j^k) and ip(k,j) is considered: it grows from at the point Wjj;^ to 
d at the point Vj or Vk, respectively. 

Next we describe how the 6 coupUngs involved look like; for simplicity we will refrain from in- 
dicating in the boundary conditions the dependence of the parameters u, Vj, w^j^k] on the distance 



The 6 interaction at the edge connecting the 7-th and k-th half line (of course, for j,k € h 
such that k e A'; only) is expressed through the conditions 



Nj =[k € m\Sjk 7t 0} U {A: e ml (3Z > OT + l)iTji i^OATai^ 0)} 

U {A; > OT -I- 1| Tjk + 0} for 7 e ot 
Nj ={k e ot| Tkj 9t 0} for 7 > ot -1- 1 

The definition of Nj has these two trivial consequences, namely 



(14) 



k€Nj ^ 7" e Nk 
j >m + l =^ Nj cm 



(15) 
(16) 



d. 
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(17) 



the 6 coupling at the endpoint of the 7-th half line (j € n) means 

(/r/0) = (Pu,k)(d) for all keNj, 



(18) 



Further relations which will help us to find the parameter dependence on d come from Taylor 
expansion. Consider first the case without any added potential, 

nmid) = ^o;«(0) + d ,)(0) + O(d^) , ^^^^ 

'P'(j,k)(d) = f'(j,k)(0) + 0(d) , j,keh. 

To take the effect of added vector potentials into account, the following lemma will prove 
useful: 

Lemma 3.2. Let us consider a line parametrized by the variable x e (0, L), L € (0, +00) u {+00}, 
and let H denote a Hamiltonian of a particle on this line interacting with a potential V, 

// = -^.y, (20) 

sufficiently regular to make H self-adjoint. We denote by if/^'' the solution of Hi// = k^if/ with the 
boundary values t/f^-'iO) = s, t//^'''iO) = t. Consider the same system with a vector potential A 
added, again sufficiently regular; the Hamiltonian is consequently given by 

d 



HA = \-&—-Aj + V. (21) 

Let 1/'^'' denote the solution ofHA^jJ - k^tfr with the same boundary values as before, i.e. 4>^a^0) = 
s, (^^'''(0) = /. Then the function ij/^J^ can be expressed as 

ifr'/ix) = s^^'^^''>^' ■if/'''(x) for all x€{0,L). 
Proof. The aim is to prove that 

-ij/'-'" + ViJ/'-' = k^i^'-' A i^'''{0) = S A i^''''iO) = t 

implies 

andi/'^''(0) = s, ifj'^'iQ) - t. The last part is obvious, since the exponential factor involved is equal 
to one, hence it suffices to prove the displayed relation. It is straightforward that the Hamiltonian 
Ha acts generally as 

Ha = — ^ + 2BA— +M' + V. 
ax^ ax 

We substitute e^-t ^^^^^ ■ ij/'-' for 1^, obtaining 

[Ha (e«r^«'^ • r'')] (X) = (e«r^fe)d^ . ^^.') (x)+ 

+ 2Mix) ^ (e^r'*fe)'k . ^s,t^ + {&A'ix) + Aixf + Vix)) e'^/o''*fe)'k . ^s,t^^^ 
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Now we express the derivatives applying the formula A(z)dz - A{x). Most of the terms 
then cancel, it remains only 

[Ha (e'^r^(''>'»^ ■ if,'-')] (x) = e«r^(''>'»^ ■ {-^•'"(x) + V(x) ■ r''(x)) . 

Due to the assumption -i//^''" + Vi//^'' = k^if/^'' , we have 

[Ha (e«r^fe)d^ . ^^.')] (x) = ;t2e«r^fe)i^ . ^^.'(x) , 

what we have set out to prove. □ 

The lemma says that adding a vector potential on an edge of a quantum graph has a very 
simple effect of changing the phase of the wave function by the value ^ A(z)dz. We will work in 
this paper with the special case of constant vector potentials on the connecting segments of the 
lengths 2d, hence the phase shift will be given here as a product of the value A and the length in 
question. 

Lemma [l!2]has the following very useful consequence. 



Corollary 3.3. Consider a quantum graph vertex with n outgoing edges indexed by I, ... ,n and 
parametrized by x E (0, Lj). Suppose that there is a 6 coupling with the parameter a at the vertex, 
and moreover, that there is a constant vector potential Aj on the j-th edge for all j € h. Let i/zj 
denote the wave function component on the j-the edge. Then the boundary conditions acquire 
the form 

,j/j(0) = i/r^(0) =: iff(Q) for all fkeh, (22) 



m, (23) 



where i^/O), tff'jiO), etc., stand for the one-sided (right) limits at x — 0. 

In other words the effect of the vector potentials on the boundary conditions corresponding to a 
"pure" 6 coupling is the following: 

• the continuity is not affected, 

• the coupling parameter is changed from a to a + B ^j- 

Proof. Consider first the situation without any vector potentials. If tf/j, j e «, denote the wave 
function components corresponding to this case, the boundary conditions expressing the S cou- 
pling have the form i.e. 

i/,0(0) = i/rO(0) -. i//°(0) for all fkeh, 

, (24) 
S^LifA" (0) = aAO). 

If there are vector potentials on the edges, Aj on the j-th edge, one has in view of the previous 
lemma, iA;(x) - e^^i'^t//\x), i.e. 

^5(x) = e-«V^;(x), 

<'(x) = ^ (e-«^->,(x)) = e-«-^^-iA/(x) - RAj ■ ^-'^'^^Pjix) . 
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Thence we express tf/jiO) and i/f^ (0): they are equal to 



«A°(0) = .A>(0), 
iA°'(0) = ^/(0)-6A,-,A,(0); 

substituting them to ( |24] i we obtain 

ilfjiO) = i/,4(0) 1/^(0) for all j,keh, 

n 

^(-A-(0)-BA,-^/0)) = a.A(0). 

./=i 

The first line expresses the continuity of the wavefunction in the vertex supporting the 6 coupling 
in the same way as in the absence of vector potentials, whereas the second line shows how the 
condition for the sum of the derivatives is changed. With the continuity in mind, we may replace 
i/'j(0) by (/'(O) obtaining 



7=1 



7=1 . 



which finishes the proof. □ 

Recall that approximating we are constructing supposes that constant vector potentials are 
added on the joining edges. If an edge of the length 2ii joins the endpoints of the j-th and ^-th 
half line, there is a constant vector potential of the value A(j_k){d) on the part of the length d closer 
to the j-th half line and a constant vector potential of the value A^i! j-,(d) = -A^jj^-^id) on the part 
of the length d closer to the fc-th half line. With regard to Lemma the impact of the added 
potentials consists in phase shifts by d ■ A(j^k){d) and d ■ A(kj)id). Let us include this effect into 
the corresponding equations, i.e. into ^T9\ : 

^a,k)(d) = e«''^"«(^|^,i)(0) + «'^[m)(0)) + 0(d^) , 

cp[j,^(d) = e«'''^<'->; .,,(0) + 0(d) , j,keh. ^^^^ 

The system of equations ( [T7| i, ([TSj, and ( |25] l describes the relations between values of wave 
functions and their derivatives at all the vertices. Next we will eliminate the terms with the 
"auxiliary" functions (p^j^^ and express the relations between 2n terms t//j(Q), if/'jiQ), j e n. 

We begin with the first one of the relations p5| ) together with the continuity requirement ([T8|, 
which yields 

d v[.,^(0) = e-«''^«"f /O) - ^|;,„(0) + 0(d^) . (26) 

The same relation holds with j replaced by k, summing them together and using the second of 
the relations (TT) we get 

{2 + d wy^) Viikm = e-«''^""^/0) + e-«''^<'-'.At(0) + O(d^) . 



We express (py it]{Q) from here and substitute into the first of the equations p5| ); using at the same 
time the first relation of (fTSl we get 
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and considering the second of the equations p5| ), we have 

if/j(0) + e'^"'('^o*-'4«..))^/,^(0) + 0(d^) 



2 + d- W{j^k] 



Since the values of vector potentials are supposed to have the "antisymmetry" property, A(^k j){d) = 
-Ay.k){d), we may simplify the above equation to 



^,(0) + e2B''Ai*^,(0)+ 0(^/2) 



2 + d- Wy^k] 

Summing the last equation over k e Nj yields 



(27) 



#Ni ■ iA,(0) = iA/(0) ■ > z ; + > i: , ^ +d-} ■ 



(#A^y denotes the cardinality of Nj), and with the help of the second of the relations ( [T8[ l we arrive 
at the final expression, 



dv 



— ; "A/CO) - > — ; — ; 



+ Oid^) . (28) 



Our objective is to choose Vj(d), wy^k](d) and A(j^k){d) in such a way that in the limit c/ ^ 
the system of relations (|28| with y e « tends to the system of n boundary conditions Q. The 
lines of (|7]l are of two types, let us recall: 



k=\ 

m 



l=m+l 



j 6 m (29) 
7 = m + 1, . . .,n. (30) 



We point out here with reference to ( [l4j i that these relations may be written also with the sum- 
mation indices running through the restricted sets, namely 



l£Nj\m 



= - 2 Tkjil^k(0) + tA>(0) 



jem (31) 
y = m + 1 , . . . , n , (32) 



keNi 



since for any pair j e m, I e {m + I, ■ ■ ■ ,n} the implication Ty/ ?t => / e A^y holds, see also 
Eqs. ([T5]l, ([^. 
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When looking for a suitable dependence of Vj(d), wyi:](d) and A^i^^id) on d, we start with 
Eq. ( |28] l in the case when j >m+\. Our aim is to find conditions under which ( |28| ) tends to (|32]i 
as — > 0. It is obvious that the sufficient conditions are 



Urn 



dvj + #Nj - y ^ 



keNj 

lim ^ 

d^o 2 + d ■ wiyjtj 

e"°'''*(y.t) 



e M\{0), 



dvj+#Nj - ZheN, 



eM\{0) '^keNj, 



(33) 
(34) 

(35) 



Now we proceed to the case j e m. Our approach is based on substitution of ( [28) into the 
left-hand side of ( [3T] i and a subsequent comparison of the right-hand sides. The substitution is 
straightforward, 



leNj\m 



leNj\ih 



VI + 



#Nj 1 
d 

#Ni 1 



d d 



-A 



lit) 



'At(O) 



f-i 2 + d- 



(d) 



'A/(0) • 



1 ^ e 

keN, 



2B"'A/.A)(^^(0) 



(a') 



(36) 



then we apply two identities, which can be easily proven, namely 

e2«'''^o.A)(/r^(0) _ y e2fi'''*(./«i/ri(0) ^ y e?-^^'-'<^\l/i(Q) 
r-d 2 + d ■ wijk) , 4t^ , 2 + d ■ wijk\ , 4Ti~ 2 + d ■ wisn 



leNj\m keNi 



^2RdA,i 



'MO) 

2 + d- W|,,^) 



leNjVn 



2 -H «i ■ W|/j) 



^/o)+ Z. Z ^^'^ 



keNiniii \leNic\m 



MO) 
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and obtain 



^;.(0)+ Yj Tjrtlr'iiO) 



leNi\m 



#Nj I ^ 1 1 y e' 



heN, 



l£Nj\m 



2 + d- W(^',(t) 



d 2 + d- W|;,,) 



l£Nt\m 

#Ni 
id) 



X'l 



dZ-i 2 + d- 



heN, ' 



<A/(0) 



(37) 



As we have announced above, the goal is to determine terms Vj{d), w\j^k]{d) and A(j^k){d) such 
that if — > 0, the right-hand side of ( |37j i tends to the eight-hand side of ( (3T| ) for all j e m. We 
observe that this will be the case provided 



' d d f-i 2 + d-w^ih) df-f ^'2 + 

heN: > leN, 



d 2 + d - wy^iij d 



leNj 
2&dAa.t) 



'jl 



^2SdA,. 



d 2 + d - wyjj 



+ T, 



lim 



leNt\m 

~d 
1 



2 + d- w^ik) 



^ S jk "^keNjnm, 



--y — 

d Z-i 2 + d- 



d^o 2 + d - w\ j^k] 



heN, 
€ 



= VZ e Nj\m , 



Mk e Ni 



1 



V/t e A^;- nm, leNkH Nj\m . 



(38) 
(39) 

(40) 

(41) 
(42) 



lim 

d^o 2 + d - W|/,i| 

It is easy to see that the set of equations ( |40] i for j e m, I e Nj\m is equivalent to the set ( |35| l for 
j > m + 1, k € Nj. Similarly, Eq. ( |42| ) for j e m, k € Nj n m, I e Nt r\ Nj\m is a weaker set of 
conditions than with j >m+\,k e Nj. Finally, Eq. ( |4T] l reduces for e Nj\m to p4l ), thus 
it suffices to consider ( |4T] i with k e Nj n m. 

The procedure of determination of v;(£/), ^{^^^((t/) and A(jk){d) will proceed in three steps, 
at the end we add the fourth step involving the verification of the limit conditions ( (33] l, ([34]l, 
and ( |4T] i restricted to A: e A^^- n m. 



Step I. We use Eq. ( |40| l to find an expression for W|^ /)(c/) and A^jj^id) when y e m and Z e Nj\m. 
We begin with rearranging Eq. (l35| into the form 



1 



2 + d- w^jJ^ 



dvi +#Ni-y - — ^ 

, « 2 + ■ 

heN, 



Wilh] j 



Ml e Nj\m . 



(43) 
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Since all the terms except e ^'^'''^ij') and T ji are real, we can obtain immediately a condition for 
A(jfl: We put 

^2B.A,„_| T,I\\TA\ if Rer,>o, 
1 -Tj,l\\Tj,\\ if Rer,v<0; 

it is easy to see that such a choice ensures that the expression ■ Tji is always real. The 

vector potential strength may be then chosen as follows, 



2d 



ai-e Ti, 



if Re Tji > 



^{srg Tji-n) if Rer,7 <0 



(44) 



for all j e m, I € Nj\m. We remark that this choice is obviously not the only one possible. 
Note that in this situation, namely if j e m and / € Nj\m, the potentials do not depend on d. 
Taking (|44]) into account, Eq. ( |43| l simplifies to 



2 + d- 



= (Tji) 



heN, 



Vl>m + l,j e Ni ; 



(45) 



note that jemAle Nj\m o / > m + 1 A 7 e A^;. The symbol (■) here has the following meaning: 
if c e C, then 

' |c| if Re c > , 
-Id if Re c < . 



Summing ( |45| l over /' € Ni we get 

y = YiTii) ■ +#Ni-y — 



I.e. 



f 



heN, 



hm j jeN, ^ ^ " ^^J''i jeN, 



We have to distinguish here two situations: 
(i) If 1 + YjheN,{Thi) + 0, one obtains 



■(dvi+#Ni), 



and the substitution of the left-hand side into the right-hand side of ( |45| l leads to the formula for 
1/(2 + d ■ wyji), namely 



1 



2 + d- w^jJ} 



- (Tj,) 



dv, + #N, 

1 + TjheNiiThl) 



e m, I e Nj\m . 



We observe that the sum in the denominator may be computed over the whole set m as well, since 
/1 ^ N/ ^ T/ji = 0, which slightly simplifies the formula. 



1 



2 + d- W|j;„ 



= {Til) ■ 



dv, + #Ni 
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'^j e m, I e Nj\m . 



From here one can easily express wy^i^, if v; is known. However, it turns out that vi(d), l> m + I 
can be chosen almost arbitrarily, the only requirements are to keep the expression dv; + #Ni 
nonzero and to satisfy ( |33) l, ( |34| i and ( |4T] i. The simplest choice possible is to define v; by the 
expression 

dvi + #Ni _ ^ 
1 + ZtiiThi) ~ ' 

which simplifies the expressions for other parameters. Here we obtain already expressions for v; 
and w\jj\ if I > m + 1, viz 



vi(d) 



l-#N, + Zti<Thi} 



VZ > OT + 1 . 



(46) 



-2 + 



(Tji) 



V j e m, I e A^;\m . 



(ii) If 1 + YjheNi(Thi} - 0, it holds necessarily dvi + #Ni - 0, and consequently, v/ 



(47) 
Note 



that this equation may be obtained from Eq. ([46|) by putting formally 1 + Yj'i'=i{Thi) - 0. It is 



easy to check that wijj\ given by Eq. ( |47] i satisfies (|43| in the case 1 + YjheNiiThi) = as well. 
Summing these facts up, we conclude that Eqs. ( |46] l, ( |47] i hold universally regardless whether 
1 + Yjh£N,{Thi) equals zero or not. 

We would like to stress that the freedom in the choice of vi{d) is a consequence of the fact 
mentioned in Remark 2.2 namely that the number of parameters of a vertex coupling decreases 
with the decreasing value of rank(B). 



Step II. Equation ( [39| l together with the results of Step I will be used to determine W{j^k](d) ™d 
A(j^k){d) in the case when j e m and k e Nj n m. From ( |39l ) we have 



2 + d- wij^k] 



2 + d- wii^k] 



the pairs {l,k) appearing in the sum are of the type examined in Step I, i.e. k e m, I e Nt\in)- 
Thus one may substitute from ( |46] l and ( |47] i to obtain 



2 + d- wij^k] 



^d-S 



jk 



2 TjiTkt- 



We observe that the summation index may run through the whole set h\m, because I > m+l Al i 
Nk Till = 0. This allows one to obtain a more elegant formula. In a similar way as above, we 
find the expression for Ay^k), 



Au,k)(d) = ^ arg 



d-Sjk+ Yj TjiTki 

l=m+\ 



for Re 



d - S jk+ ^ TjiTki 

l=m+ 1 



> 



(48a) 



and 



arg 



d - S jk+ ^ TjiTki 

l=m+ 1 





/ 


- 7T 

J 


for Re 

V 



d-Sjk + Yj TjiTki 

l=m+ 1 



<0, (48b) 
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and for w^j^k]. 



1 



2+d- Wy^k} 



d ■ S jk+ ^ TjiTki 

l=m+ 1 



(49) 



Step HI. Substitution of the results of Steps I and II into Eq. ( [38| l provides an expression for Vj(d) 
in the case when j e m. A simple calculation gives 



Vj{d) = Sjj - ^ - X (Sjx + ^ Z T^fl^ + ^ Z (1 + ^Tji)){Tj,) . (50) 

^-l \ /-m-l-l ' 



#A0 

d 



k=\ ' l=m+l ' /=m+l 

Since 5 is a self-adjoint matrix, the term S jj is real, thus the whole right-hand side is a real 
expression. 

Step IV. Finally, we verify conditions ( |33| l, ( (34] i, and ( |4T] ), the last one being restricted to e 
Nj n m. This step consists in trivial substitutions: 



(B3b : Hm 



dvj + #A^,- - y ^ 



lim 1 = 1 € M\{0) V /■ > m + 1 



34b : Um 



1 



(pTli : Hm 



d->0 2+ d ■ W^j^k] d^O 

1 



= Hm<r,y) = inj) e M\{0) V; > m + 1, ^ e A?,- , 



d^o 2 + d ■ W{j^k} ''^0 \ 
Sjem, k&Njf\m 



iim/«f • 5^, + ^ r,7r„ = Yu ^J'^") ^ 



l=m+ 1 



/=»!+ 1 



4. The norm-resolvent convergence 

In the previous section we have shown that any vertex coupling in the center point of a 
star graph may be regarded as a limit of a certain family of graphs supporting nothing but 6 
couplings, 6 interactions and constant vector potentials. The parameter values of all the (J's and 
vector potentials have been derived using a method devised originally in ||3]|3 for the case of a 
generalized point interaction on the line. The aim of this section is to give a clear meaning to 
this convergence. Specifically, we are going to show that the Hamiltonian of the approximating 
system converges to the Hamiltonian of the approximated system in the norm-resolvent sense, 
with the natural consequences for the convergence of eigenvalues, eigenfunctions, etc. 

We denote the Hamiltonian of the star graph F with the coupling (j6| at the vertex as H'^^ 
(referring to the approximatet/ system), and H^^ will stand for the approximating family of 
graphs that has been constructed in the previous section. Symbols R^^(k^) and R^^(k^) will then 
denote the resolvents of H^'^ and H^^ at the points from the resolvent set. Needless to say, 
the operators act on different spaces: R^'^{k^) on L?(G), where G = (R^)" corresponds to the star 
graph F, and R^^(k^) on L^iGd), where 



Gd = (M+)"e(0,£/)^''=' 



(51) 



Our goal is to compare these resolvents. In order to do that, we need to identify R^'^(k^) with the 
orthogonal sum 

R'^\k^) = R^\k^)®0, (52) 
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where is a zero operator acting on the space {(0, d)^!-' which is removed in the limit. 

Then both the operators R^'^(k^) and R^^(k^) are defined as acting on functions from the set 
which are vector functions with n + components; we will index the components by the 

set 

I ^h\J{(l,h)\len,heNi} . (53) 
In this setting we are able to state now the main theorem of this section and the whole paper. 

Theorem 4.1. Let Vj, j e h, , wu^k] j ^ h,k e Nj and A^^''^\d) depend on d according to ( |50| l, 
( |46] l, ( |49| l, ( |47| i, (|48jl and respectively. Then the family H'^^d) converges to H^'^ in the 
norm-resolvent sense as d — > 0+. 

Proof. We have to compare the resolvents R^^{k^) and R^^(k^). It is obviously sufficient to check 
the convergence in the Hilbert-Schmidt norm, 



Rj°(k^) - R^'^ikHl -^0+ as d- 



0+ 



in other words, to show that the difference of the corresponding resolvent kernels denoted as 
^^^''^ and Q'^'^''^, respectively, tends to zero in L^{Gd © Gd). Recall that these kernels, or Green's 
functions, are in our case matrix functions with (n + YIj^i ^j) ^ (" + ZijLi ■^;) entries. We will 
index the entries by pairs of indices taken from the set I (cf. (|53|l). 

The proof is divided into three parts. In the first and the second part we will derive the 
resolvent kernels 0^°''' and Q^^'"^ , respectively, in the last part we compare them and demonstrate 
the norm-resolvent convergence. 

/. Resolvent of the approximated Hamiltonian 

Let us construct first for the star-graph Hamiltonian with the condition (|l| at the vertex. We 
begin with n independent halflines with Dirichlet condition at its endpoints; Green's function for 
each of them is well-known to be 



Qi,K{x,y) = 



sinh Kx^ e 



where x< :— min{jc,y), x^ :- max{x,y}, and we put S>k - k assuming conventionally Re/c > 0. 
The sought Green's function is then given by Krein's formula 1 16 App. A], 



R^\k^)^R"'{k^) + J]^J>(k'){^i{l^r) 

j,i=i 



L\(S.*)") 



(54) 



where R^'(k^) acts on each half line as an integral operator with the kernel ^b^^, and for <pj{k^) 
one can choose any elements of the deficiency subspaces of the largest common restriction; we 
will work with {<p j{k^){x)^ ^ - Sjis^"^' where the symbol x stands here for the vector {x\,. . . , x„) e 
(M+)". Then we have 



R^'^ik^) 



( «Ai 



4fn 



Xi 



r+oo 

Jo &e,K(xi,yi)tJ/i(yi)dyi 

r+oo 

Jo &aAXn,yn)l^n(yn)dyn 
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+ ;^^,.,(^')(e-^\<A/(3'/)),,(, 







which should be satisfied for any (iffi, . . . , tf/„)^ e ©"^j L^(K^). We observe that for all j e n, 
the j-th component on the right hand side depends only on the variable xj. That is why one can 
consider the components as functions of one variable; we denote them as gjixj), j e n, in other 
words, 

( ipi \( xi ^ 



The functions gj are therefore given expUcitly by 

Jo Jo 



Since the resolvent maps the whole Hilbert space into the domain of the operator, these functions 
have to satisfy the boundary conditions at the vertex. 



where 



AjhghiO) + Bjhg'hiO) = for all jen. 



(55) 



h=l 



/i=l 



A = 



S 

j(n-m) 



-B = 







Using the explicit form of @t,K{xh,y) and the equality ^^g^^'''^ 



Xk=0 



= e , we find 



n p+co 

gm^yAuik^) e-'yifri(y)dy 
tt Jo 



and 



e-'>;/fH(y)dy - kV AM(k^) e-'y^i(y)dy. 
jrf Jo 

Substituting from these two relations into (|55| we get a system of equations. 



n p+ca / K « ^ 

2 j Y.'^jhAhiik^) + Bj, - Kj^BjhAhiik^) 
1=1 Jo ,,^1 



e-'yi/„(y)dy - 



with j e n. We require that the left-hand side vanishes for any il/\,iJ/2, ■ ■ ■ , ^n', this yields the 
condition AA + B- kBA - 0. From here it is easy to find the matrix A{k^): we have (A - kB)A - 
-B, and therefore 

A(k^) = (A-KBy\-B). 
Substituting the explicit forms of A and -B into the expression for A, we obtain 



A(k^) 



S + kI^'"^ kT 



j{n-m) 



j(m) J 





or explicitly 



A{e) 



{S + kI^"''> + kTT*)'^ {S + kI^""^ + kTT*)-^T 
T*(S + kI^"'^ + kTT*)-^ T*(S + kI^'"> + kTT*)-^T 
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provided that (S + kI^'"^ + kTT*)-^ is well defined. To check that the matrix S + + kTT* is 
regular, we notice that 

S + + kTT = 5 + Kif""^ + TT) , (56) 

where the matrix + TT* is positive definite and thus regular, and the value k may be chosen 
arbitrarily with the only restriction Re a- > 0. Consequently, it suffices to choose Re k big enough 
to make the matrix K{f'"^+TT*) dominate over S , which ensures the regularity of S +k(I^"'^+TT*). 

Having found the coefficients Aji{k^), we have fully determined the Green's function of 
the approximated system. Recall that it is an n x « matrix-valued function the (j, l)-th element of 
which is given by 

&^'/^^y) = ^ /^"^^-^<"'"' + Aj,(k^) e-"e- ; (57) 

we use the convention that x is from the j-th halfline and y from the Z-th one. The kernel of the 
operator R^'^{k^) is according to ( |52] i given simply by 











. 



(58) 



i.e. all entries of Q^'^^'^ except for those indexed by j, I e h vanish. 
//. Resolvents of the approximating family of Hamiltonians 

Next we will pass to resolvent construction for the approximating family of operators H^^. As a 
starting point we consider n independent halflines and Y/'j^i lines of the length d with constant 
vector potentials Ay i)(d), both halflines and lines of the finite length are supposed to have Dirich- 
let endpoints. We know that the Green's function is 0iiKix,y) = /f ' sinh A-ji[:<e in the case of 
the halflines. The Green's function in the case of the lines of the length d will be found in two 
steps. We begin with a line without vector potential and with Dirichlet endpoints; the Green's 
function can be easily derived being equal to 

sinh Kx^ sinh K{d - x>) 

§&K(x,y) = r-^— 1 ■ 

K smh Kd 

The Hamiltonian of a free particle on a line segment acts as if a vector potential A is added 



it changes to (-6;^ - . Using Lemma 



3.3 



it is easy to check that 



d .\' d2 



where U is the unitary operator acting as 



m)(x) ^ e'^^'i/^ix) . 



If we denote Hq - --^ and Ha = (-Kg^ ~ ^) ^ we see that 

(Ha - A)-' = (UHoW - Ay' = (U(Ho - A) Wy' ^ U (Hq - Ay' W , 
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so the corresponding resolvents are related by the relation analogous to (|59|. This yields 



' e'^^'- si"h.-^<sinh4J-x>) ^-My^^^^^ ^ 
K sinh Kd 



thus the sought integral kernel is equal to 



rA ,A _ „BA.v Sinh/fx< sinh;^(^-^>) 
K smh Kd 

Now we can proceed to the derivation of the complete resolvent R^°{k^) which will be done 
again by means of the Krein's formula. The situation here is more complicated than in the case 
of the approximated system; recall that R^°{k^), as well as H^'^{k^), acts on the larger Hilbert 
space L^(Gd), where Gd has been defined in ( |5T| i. Moreover, the application of Krein's formula 
means that we have to connect all the line segments using the appropriate boundary conditions, 
i.e. we must change boundary conditions at n + 2 'Z/j=i Nj endpoints, specifically n belonging to 
n half lines and 2 Yj'j=\ belonging to segments of the length d. Thus the index set for 

the indices in the sum on the right hand side of the formula has n + 2 Y/j=i Nj elements; we will 
index them by the set 

J = n U [(l,hf \ lehjie Ni] u[(l,h)''\ leh,heNi] . 

The elements of n correspond to changed boundary conditions at the endpoints of the half Unes, 
and the elements of the type (/, /i)" and (/, h)'^ (h e Ni) correspond to changed boundary conditions 
at the endpoints of the segments of the length d which are connected to the endpoint of the /-th 
half line. If we denote by the symbol R^'^ik^) the resolvent of the system of the n + Nj 
decomposed edges with Dirichlet boundary conditions at the endpoints, Krein's formula for this 
pair of operators has the form 

J,L£l 

The role of the superscript d in the lambda symbols is to distinguish them from Aji that have been 
used in Eq. ( |54] l for the resolvent of the approximated system. The functions 0^ (J e I) may be 
chosen, as before in the case of the approximated system, as any elements of the corresponding 
deficiency subspaces of the largest common restriction. Note that each function (/>^ has n+2"=i Nj 
components indexed by elements of the set J = n U {{l,h)\ I eh, he Ni]. It turns out that a 
suitable choice is 

5,xe— i for jeh,LeI, 

= efi^<"rt'-'>5(,_;,)£sinh^X(,,;,) for leh,heNi,LeI (61) 

I'P'lj.yik^X^)^ = efi^""'^<'.")5(, sinh K(d - ^(,,,)) for I e h, h e N, ,L e I , 

where the symbol x denotes the vector from Gd with the components indexed by I. We remark 
that if 7 e n, 0^ is independent of d and equal to the corresponding function chosen above in the 
case of the approximated system. 
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If we apply the operator (|60| to an arbitrary e L^{Gd), we obtain a vector function 

with n + 'T!]=\ Nj components indexed by I, we denote them by gj (j € n) and g(/ /,) with I e h,h e 
Ni. As in the case of the approximated system, a component depends on xy only, thus each gj 
can be considered as a function of a single variable. A calculation leads to the following explicit 
expressions for gj,j e « and g(ij,),l e fi,h e N/; for better clarity we distinguish the integral 
variables on R'*' and on (0, d) by a tilde, i.e. y e M^, y e (0, d). 

X+oo /-*+co 
g,Axj,y)il^j(y)dy + Y^Ajj,(k^) c'^y ■ ^j,(y)dy c-"""' 

+ Z Z f e-'^"'^-'^'&inhKy-i/,g,,h'){y)dy 

/'=1 h'eN,' ^ 



(62a) 



Jo 



"(j^(/,/i),5')'A(/,/o(5')d5' 



+ e 



. sinh/fjC(,,/, 



/=i 



/' = 1 /I'SW,, ^ 



e '^'^c sinh /fy ■ (/'(/', ft' )()') 



e-H''"''-^sinh<flf-y)-(A(/',/,')(j)dJ 



+ e'""'"'"'" ■sinh/f(t/-X(,,/,))- 



+ Kimi'h'Y^^^'^ e-"^^"' sinh^flf - y) ■ il>o',h')(y) dyj 



By definition the function (gj)jei belongs to the domain of the operator H^^, in particular, it has 
to satisfy the boundary conditions at the points where the edges are connected by 6 interactions 
and 6 couplings. Step by step we will write down now all these boundary conditions; this wiU 
lead to the explicit expressions for the coefficients Aj^(k^). 

Step 1. The continuity at the points W\ j^k] means 

giuhm = s^m (63) 



(62b) 
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for nil I en, he N,. Since ^'*(0,y) = for all y e (0, d), it holds 



8(i,h)(0) = sinh Kd- 



/'=1 /I'eW,, ^ 



Jo 



the expression for g(i,j){Q) is similar, just the positions of / and h are interchanged. Since Eq. ( |63| l 
must be satisfied for any choice of the function = (i/ry) yj=j, the following equalities obviously 
hold for V/ e h,h e Nr. 



^ (1,11)''/ ' ^ iKVr ' 2 



V/ e h , 

V/' eh,h' e Ni' , 



(64) 



In other words, all the coefficients /l^^ with J e I(k^) are symmetric with respect to an 

interchange of I and h. 



Step 2. The sum of derivatives in points is 
for all I e h,h e Ni. We substitute 



(65) 



dx 



sinh K(d -y) 
sinh /fc/ 



and 



(e'^'^^sinh/cjcVl = /c 



into Eq. ( [65| ) and require the equality to be satisfied for any *P - {ipj)jei- In the course of the 
calculation, the outcome of the Step 1 is also used. As a result, we find how the coeflicients 
"^m"/'^^^ ^-^ ^ -^^ expressed in terms of ^f,,,pjik^) and ■^fi^j-^ojik^)- 



U.hYf 



2k cosh Kd + w^i^i,} sinh Kd 



(66a) 



(i,h) {i',h') 2k cosh Kd + w\i,h] sinhATii 



(66b) 



2k cosh Kd + W{ih\ svahKd 



vr e h, h' e N,, (66c) 
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for the indices I e n,h e Ni. 

Step 3. The continuity at the points Vj requires 

gj(0) = g(jj,^(d) (67) 
for all eh, he Nj. Since ^8-^(0,}') = for all jc e M+ and g^^'Hcl) = for all y e (0, d), it holds 



+ Z Z pe-«^<"'"'^inh^^^(,,„ 



e-V,.(y)dye-"^ 

)0') dy 

+ A'.^,,,^,/k^) e-«^"' '"'^ sinh ^(flf - mi',h')(y) dy) - (68) 



and 



f' y-»-f-00 

/=i 

+ Z Z e-«V.„',Vsinh^5)^^(,,,,)(5))d5) 



(69) 

The relation (|67]i should be satisfied for any choice of *P = ((/r y) y^j, hence we obtain the coeffi- 



cients A'l^j^^„j(J?) with J e I expressed in terms of A.''jj(k^) in the following way 

4/.)»(/',/.')«(^') = ■ ^KVMf^^'^ e «' e ' (70b) 

for j e h,h e Nj. We also return to the result of Step 2 - we substitute there for A''^. ^^^^^ j(k^) the 
expressions that we have just obtained arriving thus at 

Ai , w ., {k^) = ^— ■ U-^''^'"'^ Ai, (k^) + e-'^"''*"'" A^ ., (k^)] 'if eh, 

J 2/c cosh /ft/ + sinh sinh/cii ^ -' -' ' 

(71a) 

.rJk^) = 1_ .(e-B''Aw4 ^, Jk^) + e-^-^^w'if,,, „Jk^)] 

(imi'.h'r > cosh /ft/ + W|,,ft) sinh /fc/ sinh/fc/ ^ '(' ■'' '"^ " 

Ml' eh,h' eN,, , (71b) 
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1 



2k cosh Kd + W(/ /,) sinh Kd sinh Kd 



-MA, 



hii'M'Y 



Ml' eh,h' eNi, (71c) 



for I eh, h e Ni. 



Step 4. In this step we examine the sum of derivatives at the points Vj, i.e. at the junctions 
of the halflines and connecting segments. Since the connecting lines support constant vector 



potentials, one has to rewrite the original condition into the form derived in Corollary 3.3 Note 



that the variable on the connecting segments is considered in the ingoing sense, thus the sign of 
the potentials Aqj,-, {h e Nj) has to be taken with the minus sign. The resulting condition is 



heN; 



heN; 



8ji0) 



(72) 



for all Jen. 

The way how to proceed in this step is essentially the same as in previous steps, only the 
calculus is slightly longer. With the aid of the formulae 



dx 



dx 



x=d 



x=0 

sinh Kd 



sinh Ky e 



-BAv 



and 



^ ^ \x=d 

{c'^^-'sir)hK{d-x))'\ ^^-Kt^^'', 



used in Eq. ( |72] i, we arrive at an expression containing t//i . . .,4r„ that should be satisfied for any 
choice of ^l' = {4'j)jei- This yields the following three groups of conditions; 



5jj. - xA'^jfik^) - [k cosh Kd + S,Ay,h) sinh/fd) ■ e'^'''*<'"'^^^.^j„^,(/t^) 



heN: 



heN, 



Ajj,(k^), (73a) 



sinh Kd 



K^%j,/k^) - X ('^cosh^flf + BAo;,,) smhKd) ■ e'''^"'"4.,,)„(,,,)„(^') 



heN; 



heN; 



heN; 



^%j,f(k^), (73b) 
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heNi 



heNi 



MM')'' 



(n- (73c) 



We use the equalities ( fTO] ) and ( |7T| l to eliminate all terms of the type ^'Iji^fjik^) and A'^^^^^yjik^), 
J e I. In this way we obtain three independent systems of equations for Ajj,{k^) e n), 
4,,,,,)o(^') ^' eh,h'e Nr) and I' eh,h'e N,,): 



,j T cosh/fc/ ^ T x— 1 



^2 AjjXk^) + e^'^''^^fWlj,(k^) 



sinhATfi 



hsN, 



sinh /fd 2/<- cosh Kd + W| jj,) sinh Kd 



-vjA'^jXk'), 
(74a) 



sinh Kd 



cosh /ft/ 



E 



j(l.h'r- ' "^MLh'f^ ' 



sinh /cd 2k cosh /ft/ + W|y,/,) sinh Kd 



vj^%j,^o(k'), (74b) 



cosh.^.. V. .2 + ^'""""'4/,///^') 



J ^ cosnra J T v 



sinh /cc/ 2/*: cosh Kd + W|^_;,) sinh /ct/ 

IT^ ■ A ;S ir^ r^6j,,dhh' = v,-4,, Xfc') • (74c) 

ih k-d i—i 7.*-cosh kd + wi v/.i sinh icd 



sinh Kd ^ 2k cosh Kd + wii u sinh /cc/ 
Let us focus, e.g., on Eq. ( |74a[ ), which can be rewritten in the form 

cosh Kd 



z 

/7=1 



/f + K#Ni 



^ sinh A-ii sinh Kd f—^ 2k cosh /c^^ + w, , r,, sinh Kd 



+ V/ 



sinh Kd 2k cosh /C(i + sinh /ct/ 



4'(^') = ^77' (75) 



for all j, f e n; the symbol ;(fAf^(/!) is equal to one \fhe Nj holds and zero otherwise. As we will 
see within a short time, it is convenient to introduce a matrix M^/ the (j, h)-th element of which is 
defined by 



[Md]jh = Sji, 



K + K#N 



cosh Kd 



— y 



^ sinh /<•£/ sinh Kd ' 2/i- cosh Kd + w, , 7,, sinh Kd 



+ V, 



-XN,{h) ■ 



sinh /ct/ 2k cosh /ft/ + W(j./,| sinh Kd 



■ (76) 
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We also rewrite the set i" as a union, 



J = « U J"" U J"'' , 

where - ((/,/!)°| / € h,h € A^/} and J'' - i^(l,hY\ I eh, he A^;}, and define the symbols 
h-xfik^) for X,Y e {n, J"''} by the relation 

e.g. K^^^{k^) = [A''.,,, .,Jk^)\ ■ Obviously, the matiix A^^'''(k^) has the block structure 













A^jAk') 




AX^ik^) 





We observe that the system of equations ( |75] l is nothing but 



and therefore h^^f\k^) = (M^)-', or in the components 



(77a) 



The matrices A^^^(k^) and A^^j(k^) can be found in a similar way. We start from Eqs. ( |74b| i 
and ( |74c[ ) and arrive at 



and 



To obtain expressions for A^^'fc^^) and A^ik^) = J^,J'') we substitute ( |77al ), (|77b]) and 



sinh /cJ 2/c cosh Kd + wii\h'} sinh Art/ 
( |77c| l into Equations ( |70] i and ( j71| i which gives 



(77b) 
(77c) 



A 



^ " ■ ^ii^h^ ill' ' 

j-B£/A(,,,,) ^ . gBi/A|,/ j/i 



^(WC/',/,')"*^^ ) ^jjjj^2 2/f cosh /ft/ + sinh Kd 



(77d) 
(77e) 

[(Mrf)-']^^, , (77f) 



(/,/!)/ 2/f cosh /ct/ + W|//,| sinh /cc/ sinh 



(77g) 
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2k cosh Kd + W|/_/,) sinh Kd sinh^ Kd 



-MAu 



+ e 



[(M.)-'],J, (77h) 



K 



1 



2k cosh Kd + w\ih\ sinh Kd sinh /cJ 



sinh /cc/ 2k cosh A-t/ + W|;'_//) sinh Kd 



h)(l',h') 



(77i) 



Once we compute the elements of (M^)"' expHcitly, we will have fully explicit formulae for 
A^°'''(k^) we need. We start from the matrix M^/ itself. We take the formula (|76|, substitute 
there the expressions for Vj(d), wij^k}{d) and AQ^ii)(d) that have been obtained heuristically in 
the previous section and apply Taylor expansions to appropriate orders. A slightly laborious 
calculation leads to the formulae written below; note that the structure of the expression for the 
(j, h)-th element of the matrix depends on whether j, h belong to m or to h\m: 



[Md]jh = Sjh 



/=m+l 



l=m+ 1 



(78a) 



[Md]jh = --Tjh + (d) for jem,h>m+l, 

[Md]jh = --Thi + (d) for j>m+l,hem, 
d 



[Md]jh =Sjh\K+ - \ + (d) for j, h>m+l. 



The matrix M^/ has thus the form 







a 


[k + i) !'•"-'"'> 



(78b) 
(78c) 
(78d) 



where (d) on the eight-hand side represents a matrix nxn the all entries of which are of order of 
(d) as — > 0. 

Our aim is to find the inverse of M^. For this purpose, we denote the first term on the right- 
hand side, the principal one, as M/j p and notice that if M^ p is regular, then 

[Md]-' = iMd,P + (d))-' = [Md,pU + [Md,p]-\d))l' = 

= [/ - [M,,,p]-'(d)] [Md,p]-' = [Md,p]-' - [M,,,p]-\d)[M,,,p]-' . (79) 

Moreover, if [Md^pT^ - (1) it obviously holds [M^]"' - [Mj,/>]"' + (d); in other words, under 
certain assumptions it suffices to find the inverse of Mj p. 
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Since the matrix p has a block structure, one can find M^i p in the same block structure. 



This in other words means that we are looking for a matrix 







A^3 


N4 



such that the relation 





N2 \ 


( S +kI^"'^ + \TT* 


1 7- ^ 


/ /(«!) 





A^3 




1 'X'* 




"I 


j(n-m) 



holds true. It turns out that 

[Md,p]-' = 



1 



(\+Kd)- 



_j_(5+^/('") + ^rrT'r 

T*(S + kI^'"^ + Y^Try^T + T^/*"""'^ 



provided the matrix S + kI^'"^ + jfj^TT* is regular. Since - k + (d), we may proceed in the 
same way as in (|79|, and we obtain 



(s + /f/*'"' + — ^rr*) ' = (5 + /f/*'"* + kTT* + (d))-^ = (5 + + kTT*)-^ + (d) , 

\ \ + Kd I 

if the matrix S +kI^'"^ + kTT* is regular. However, the regularity of this matrix has been discussed 
and proven for an appropriate k at the end of the part devoted to the approximated system, see 
Eq. ( |56] l and the paragraph following it. It follows that M^^p is regular as well, i.e. the condition 
on the regularity of Md,p in ^T9\ is satisfied. Hence 



[Md,p]-' 
and together with ( [79| l we have 

It is important to notice that 



(S + K^"^ + kTT*)-^ (S + kI'-'"^ + kTT*)-^T 
T*(S + kI'^'"> + kTT*)-^ T*(S + kI^'"^ + kTT*)-^T 



(s + ki^'"^ + KTry^T 



+ (d), 



+ (d). 



[MJ-i=(l) fort/^0+. 



(80) 



Combining the above result with Eq. ( |77a[ ), we can conclude that 



nn ^ ' 



T^iS + kI^'"^ + kTT*)-^ T*(S + kI^""' + kTT*)-^T 



+ (d). 



hence 



(81) 



Having the coefficient matrix we can determine the resolvent kernel. First we introduce 
symbol ^ = { (/, h)\ I e h,h e Ni} (i.e. I - h U J'), then we employ a notation similar to the 
case of the matrix A^^'''{k^) and its submatrices. We introduce symbols 0xYk ^'^^ ^'^^ P^^'^^ X,Y e 



{h,J} to denote the blocks = 
the structure 



JeX.LeY 



; then the integral kernel Q^^''' of R^^{k^) has 
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(82) 



for x,y e Gd- Using ( |58| ) we can write the diiference in question as 



- 0R 



' ^Ag.rf ^Ad 






^Ag,i/ 



(83) 



///. Comparison of the resolvents 

To make use of the above results we compute first explicit expressions for all the entries of 
Q^''^(x,y), up to the error term in the lambda coefficients indicated in ( |8T| ). They may be derived 
from Eqs (|62]| together with ( |6T] i: 

Qu'n'^^pyi'^ - ^n' ''"^""^ + 4(^^) e-'e--^/ for (84a) 



'■"''"''■"''■K''.A')»(^^^ sinh/fy(,V/) + sinh^t/ - y(,v/))] 

for jeh,(l',h')ej, (84b) 



B^w/^-^C.*)'^/) = e^^<w^™ ■ [4,,)„,.,(^') sinh^^(y,) + 4,).^,(^') sinh^(fl' - X(,,a))] ■ e"^^'/ 

for il,h)ej,j' eh, (84c) 



,. sinh/fx^ sinhA'ft/ — x>) „. 



Ag,rf 

&K.mil',h')y-V..n.,.,n„ -U,«A...W- /fsinh/ft/ 

+ e«^<'-«"<'-« ■ sinh ■ e-'^^ii'y>yi''y> ■ [4,,)„(,,^,)„ sinh ^yf,,,^,) + if, /„0(,,,,,). sinh ^(of - y^r,,,,^)] 



+ e 



BA,„„.„.,„ . sinh K{d - xg^h)) ■ e-^^<'' •*'»•<''•'■'' ■ [if,,;,).(,,,,)0 sinh>o'(i'.A') + if,,ft).(,,,,-). sinh^(flf - y(,,,^,))l 

for (l,h),(l',h')ej. (84d) 

Now we are able to compare the entries of ^^^'"'(x, y) given by ( (82] i and 0^'^'^(x, y) as speci- 
fied in (|58]l). We begin with the upper left submatrix « x « of ( |83| ). From the expressions for 
and ^g^^^l, cf. (|57]) and ([84a|, we have 



sinh /cx< e 



sinh /(•x< e 



= {A''jj,(k^) - Ajj,{k^)) e-'^^ie-'y'' = [h-'^fik^) - A^''(/t^)] ., e^'^^^e-'^-'''' = (^^) e-^^e" 



(85) 



where the last equaUty holds by virtue of ( [8T| l. Since such estimate is valid for all j, f e «, there 
is a constant K\ independent of j, f and d such that 



holds for all j, f e h, Xj,yj' e and any d sufficiently small. 
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(86) 



Then we proceed to the upper right submatrix of ( [83) . To find a bound for the entries of 
(|77b|) and (|77cli 



@'^°/^{x,y), cf. ( |84b[ ), we substitute values of /l^^^, ^,jo(^^) ™d /l^^^, ^lyik^) that we have obtained in 



-KXj .p-6A„/ ,,/|V|/'.,/) .p6rfA(,, v) 



sinh*fy(/-,//) 



sinh Kd sinh /cc/ 2/c cosh Kd + W|;'_/,') sinh a-c/ 
It holds [(Mrf)"'j - (1) by virtue of ([80| and obviously |e"H/'j.').V|/'.,,'i 
suffices to estimate the terms in the brackets. When d is sufficiently small, it holds 



1, thus it 

sinh /f J I ' 

< 1 . As for the denominator of the second term, 
we substitute for W(/',/,') from ( |49| l or ( |47] i, depending on whether both /', h' belong to m or not, 
and we easily obtain the estimate 



because < Jq'Ji') < d; similarly 



| sinhy(t/-V(,<j,/|) | 
sinh Kd 



1 



2k cosh Kd + W(/'_/,') sinh Art/ 



= (1). 



Summing all this up, we get 



independently of j, (I' ,h') and x,}?, thus there is a constant 7*^2 independent of d such that 



(xj,y(i',h'))\ <K2& 



(87) 



for all j e n, {V ,h') e JT, -Xy e M^, ^(/'yi') e (0,c/) and d sufficiently small. 

Similarly we proceed in the case of the left and right bottom submatrices of ( [83] ), i.e. when es- 
timating the entries of Q^°jJ(X,y) and Q^^jj{x,y). As for Q'^^jJ(X,y), we substitute for ■^^i /^-^o jXk^) 
and A'l^i^^,.,(k^) from (|77d| and ( [77g| ) into ( [84c) and obtain 



(/./il 



sinh KX(^ifi) 
sinh /(-(t/ - 



sinh /ft/ L J/> 
1 



2k cosh Kd + W|; /,) sinhAc/ sinh At/ 



Using analogous estimates as in the case of @^^''lfj, above, we obtain 

thus there is a constant Kj, independent of d such that 

\0t:LAwf)\<K3^-'''' 

for all (Z, h) e J', f e n, e (0, J), y e and d sufficiently small. 
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(88) 



Finally, we substitute from (|77eJ, (|77h| and ^ for ^;',_„)0(,,,,)o(^'), <„)0(„,)^^^'^' 

Kuhnvh'f^^^^ K,h)''(i'h'A'^^^' respectively, into Eq. ([84d| and obtain 

„\ffri Ri ,■ sinh/('x< sinh x>) „. 

»B/f,(;,/!)(/',/i')'- -^C ■'')'' "(iMi'Ji')^ KsintiKcl 



+ e 
sinh 

1 . 

sinhA-ii siDhKd 2*- cosh Art/ + W(;'_/,') sinh /cii 



sinh Kd 



sinh ,/,-)) 



gBrfAy/j,,, 



sinlW(^A') /^-BJA,, 
sinh A-t/ 



+ e 



BdA|,'.t 



2*' cosh /f J + W{ih] sinh A-ii sinh /<-(i 
sinh K(d-y(i,jj,)) 



sinh A-c/ 2a- cosh /cc/ + W|/' /,/) sinh Kd 



It obviously holds 

thus there is a constant independent of d such that 



(89) 



(90) 



for all (/, h), (I', h') e JT, xijji),y(v,h') e (0, li) and any li sufficiently small. 

With the help of ( (86] l, ([87|, ( |88| l and (|90|, we may now estimate all the entries of ( (83] l, which 
will allows us to assess the Hilbert-Schmidt norm of the resolvent difference for the operators 
H^^ and Z/^. This norm can be written explicitly as follows, 



<(^^)-/^f(fc^) =x r r 

112 j^.Jo Jo 



>1 (/',A')eJ 

(/,/!)£/ / = 1 ^" 



y(Z',A')) dxydyy-,/,') 



,,yXjC(W,y/)| dx(y,)dy^v 



Zj Z L [, \^&lm(i'h' 



,(-^(/,/i)'3'(/,A')) ^m,H)Ay(v ,h') ■ 



(l,h)£l {l',h')el 
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Now we employ the estimates derived above obtaining 



X r P|/^.Je-^e-'^/pdx,dy,,+|^ Z C fl^ie-f dx,d);y,v) 
ZZ |ii:3e-^/f dxa,)d>',.+ 2 Z 1^:4!' dxa^)d)'(, 



+ 



" p+00 pd 

+ Z Z ^2 r ^"''"''^^^dx,. r ldy(,,v) 
ldx(/,ft)- I ld);y,,v) 

J,/ = l ^ ^"^J j=l (,l',h')€l (.l,h)^IJ' = l (/,/!)€/ (/',V)€ 



= 

Hence 

- /^^''('fc^)!^ = ( for ^ 0+ , 

and consequently, the Hilbert-Schmidt norm of the difference R^^(k^) - R^'^(k^) tends to zero as 
d ^ 0+ with the explicit convergence rate. Since the HS norm dominates the operator one, it 
follows immediately 

jim||<^(^)-/?f(fe2)||=0, 

therefore the resolvent difference tends to zero in I?{Gd) as J — > 0+, which we set out to prove. 

□ 
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